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Abstract. The population of compact extragalactic sources contribute to the non-Gaussianity (NG) 
at Cosmic Microwave Background (CMB) frequencies. We study their NG using publicly available 
full-sky simulations. We introduce a parametrisation to visualise efficiently the bispectrum and we 
describe the scale and frequency dependences of the bispectrum of radio and IR point-sources (PS). 
We show that the bispectrum is well fitted by an analytical prescription. We find that the clustering of 
IR sources enhances their NG by several orders of magnitude, and that their bispectrum peaks in the 
squeezed triangles. Examining the impact of these sources on primordial non-Gaussianity estimation, 
we find that radio sources yield an important positive bias to local /nl at low frequencies but this 
bias is efficiently reduced by masking detectable sources. IR sources produce a negative bias at high 
frequencies, which is not dimmed by the masking as their clustering is dominated by faint sources. 

1. Introduction 

The CMB, the dominant signal on the sky around 100 GHz, is a powerful probe of the early universe 
which significantly contributes to the establishment of the standard model of cosmology. In the 30- 
350 GHz frequency range, other signals contribute at small angular scales and modify the statistical 
distribution (initially close to Gaussian) of the measured CMB anisotropies. Three of these signals 
are associated with extragalactic sources: 

- galaxy clusters, in which the electrons of the hot ionised gas scatter off the CMB photons leaving a 
distinct spectral signature (Sunyaev-Zeldovich effect) 

- radio loud galaxies with Active Galactic Nuclei emitting through synchrotron and free-free processes 

- dusty star-forming galaxies, where the UV emission from stars heats the dust which consequently 
reemits in the infrared domain 

We will focus on the characterisation of the distribution of these two latter populations. Radio 
sources can be considered randomly distributed on the sky and are hence modeled as a white-noise 
entirely described by the sources number counts. Their harmonic "N-point" correlation functions 
(power spectrum, bispectrum, etc) are then flat and related to the corresponding moment of the 1- 
point distribution. By contrast, IR sources are highly clustered therefore their harmonic correlation 
functions require the inclusion of the spatial distribution of the sources in addition to their number 
counts. 

The primary NG in the CMB can break the degeneracy between the primordial processes 
generating the cosmological perturbations e.g. the inflations models. While most models predict 
similar power spectra, they may be distinguished at the 3-point or higher level. As an example, 
standard single-field models predict undetectably small NG while multifield models, non-canonical 
kinetic terms or vacuum initial condition etc produce potentially detectable NG. As of today, no 
definite deviation from Gaussianity has been found in the CMB data, even if there is a hint of 'local'- 
type NG in WMAP7 data [1]. The local- type NG, produced generically by several inflation models, 



takes the form : 

$(x) = <f G (x) + / NL ($ G (x) 2 - ($ G (x) 2 )) (1) 

where $ is the Bardeen potential and the subscript G denotes the Gaussian linear part. 

To study the distribution of radio and IR galaxies we use publicly available full-sky simulated 
maps of these sources [2] between 30 and 350 GHz at high angular resolution. The results presented 
in this communication are presented in greater details in [3] 



2. Bispectrum 

A gaussian field on the sphere is entirely characterised by its mean and 2-point correlation function (or 
its power spectrum). The 3-point correlation function is thus the lowest order, and most prominent, 
indicator of NG. The bispectrum is this 3-point correlation function in harmonic space : 
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A full calculation of the bispectrum at WMAP or Planck resolution (£ mSLX ~ 700 and £ ma: _ 
respectively) is computationally very intensive so one usually resort to binning multipoles. 
The bispectrum is invariant under permutation of (^1,^2^3)) it is a function of the triangle shape 
only. Three triangles shapes are of particular importance : squeezed (one side much smaller than 
the two others), equilateral and folded (flat isosceles). To plot a bispectrum in a efficient way, e.g. 
without redundancy of information, we introduced a parametrisation accounting for this permutation 
invariance, and we defined three parameters (P,F,S) based on the elementary symmetric polynomials 
o\ = i\ + £2 + £3, o"2 = £±£2 + £i£3 + ^2-^3 and 03 = ^i^2^3- The locations of the different triangles of 
constant perimeter P in the (F,S) plan is shown in Fig.l (left panel). One can then plot a bispectrum 
by making slices of perimeters and colour-coding the value of the bispectrum. 
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Figure 1. Left 
panel: Positions 
of the triangle of 
a given perimeter 
in the (F,S) plan 
of the parametri- 
sation. 

Right panel: 

IR bispectrum at 
148 GHz in our 
parametrisation 



The local- type primordial NG defined in Eq.l leaves a characteristic imprint on the CMB, which 
can be detected at the 3-point level. The bispectrum thus generated peaks in the squeezed triangles, 
and because it has a separable form [4] it has been shown that a fast estimator for its amplitude 
/nl can be built, which does not need the -prohibitive- computation of the full bispectrum. This, 
commonly-called KSW, estimator is optimal in the sense that it minimizes the \ 2 °f the fit of the 
observed bispectrum to the local bispectrum : 
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where N^ 2 ^ = ( 2 ^ 1+1 )( 2 ^|+ 1 )( 2 ^+i) ^ ^2 ^3^ - g ^ e ^^gj. Q f guch configurations on the sky 
3. Point-sources 

Radio sources are randomly distributed on the sky, so they produce a constant bispectrum b ps oc 
J S 3 ^dS, but as higher frequencies are surveyed e.g. with Planck, the IR population becomes of 
growing importance and one must account for the superposition of these two signals. From a method 
proposed in [5] , we have developped a full-sky analytical prescription to compute the bispectrum and 
higher order moments based on number counts and power spectrum. The prescription can handle 
individual populations and their superposition [3]. Analytically, at the bispectrum level, it yields : 

bf t » = aJcyt-Cf-CP- with a= J dS * (5) 

We used publicly available simulations of these two populations between 30 and 350 GHz by [2] 
to compare the measured bispectrum with our prescription. Fig. 2 shows the measured and predicted 
bispectra at 350 GHz in the case of superposition of IR and radio sources. 
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Figure 2. IR+radio bispectrum at 
350 GHz in some configurations. Black 
line is the measured bispectrum from 
the simulations, red line is the result 
of prescription, blue dashed line is the 
case where radio and IR populations are 
assumed correlated with each other. 



We found a good agreement, within cosmic variance error bars, between the prescription and the 
measured bispectrum at all frequencies and in all cases: radio sources alone, IR alone, radio and IR 
together. It is also noticeable that the IR clustering increases NG by several order of magnitudes at 
large angular scales - as compared to unclustered case. While the radio bispectrum is flat, the IR one 
exhibits a configuration dependence due to clustering (see at 148 GHz visualisation in Fig.l, right 
panel) . The IR bispectrum peaks in the squeezed triangles -as predicted by our prescription- with a 
slight flattening at high multipoles due to shot-noise, especially at the higher frequencies. 



4. Estimating the contamination to primordial NG 

We used the KSW estimator described previously (Eq.3, [4]) on the simulated maps, to estimate the 
bias of the primordial local NG parameter /nl due to foreground signals. We computed this bias 
for two angular resolution (WMAP-like, Planck-like), with and without masking sources above the 
Planck ERCSC catalog [6] flux cuts. Results are presented in Table 1&2 for radio and IR sources 
respectively. 

The relative error bars of these biases are of 1.5-2.5% for radio sources and 3 to 7% for IR sources. 
We find that for both populations the bias increases at higher maximum multipoles, as the CMB 



Table 1. Bias from radio sources to local NG parameter /nl, depending on frequency, angular 
resolution, and masking or not detectable sources 



without flux cut 



v (GH/ 


30 


90 148 219 277 350 


/ — 

/ — 
/ — 


50 -4.2 - 
700 3850 
2048 177000 


•0.0025-0.00037-0.00021-0.00027-0.00068 
2.5 0.38 0.21 0.27 0.65 
117 18 9.7 12 30 


with flux cut 


/ — 

/ — 


700 108 
2048 4930 


0.17 0.0071 0.0031 0.0035 0.0064 
7.5 0.31 0.14 0.16 0.29 



Table 2. Bias from IR sources to local NG parameter /nl 



without flux cut 


v (GHz) 


30 


90 148 


219 277 


350 


Anax — 50 
Anax = 700 
Anax = 2048 


-3.5-10" 
-1.3-10' 
-1.8-10- 


- 8 -5.9-10" t, -9.2-10- 5 - 
' 6 -0.00019 -0.0033 
~ 5 -0.0026 -0.039 


-0.0027-0.023-1.0 
-0.063 -0.55-9.0 
-0.68 -4.8 -67 


with flux cut 


*max = 700 
Anax = 2048 


-1.3-10 
-1.8-10- 


6 -0.00019 -0.0033 
' 5 -0.0026 -0.039 


-0.078 -0.74 
-0.67 -6.3 


-11 
-66 



signal plummets due to Silk damping. The radio bias is positive and maximal at low frequencies, 
while the IR bias is negative and peaks at the higher frequencies. At a WMAP-like resolution, except 
at 30 GHz, these biases can be neglected compared tho the expected error bars on /nl due to cosmic 
variance. However at Planck-like resolution, these biases become non-negligible ; while the radio bias 
is very efficiently reduced by masking detectable sources, it is not the case for IR sources where the 
bias is unaffected. Indeed the IR population, and especially the clustered galaxies, is dominated by 
faint sources much below the detection limit. 
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